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Abstract. Bound and resonance states of quantum dots play a significant role in 
photo-absorption processes. In this work, we analyze a cylindrical quantum dot, 
its spectrum and, in particular, the behaviour of the lowest resonance state when 
a magnetic field is applied along the symmetry axis of the cylinder. To obtain the 
energy and width of the resonance we use the complex rotation method. As it is 
expected the structure of the spectrum is strongly influenced by the Landau levels 
associated to the magnetic field. We show how this structure affects the behaviour of 
the resonance state and that the binding of the resonance has a clear interpretation 
in terms of the Landau levels and the probability of localization of the resonance 
state. The localization probability and the fidelity of the lowest energy state allows to 
identify two different physical regimes, a large field-small quantum dot radius regime 
and a small field-large quantum dot radius, where the binding of the resonance is 
dominated by the field strength or the potential well, respectively. 
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1. Introduction 

The availability of quantum devices with characteristic length of only a few nanometer 
allows the implementation of a number of experimental setups that put under test the 
very foundations of Quantum Mechanics. The fact that the radius of a Landau level 
is of only a few nanometers for magnetic field strengths of around 10 Tesla allows 
the verification of the Aharonov-Bohm effect in quantum rings [1], which leads to the 
presence of persistent currents [2]. This persistent currents have been measured even 
for one electron states [3]. The quantum rings are formed by only one semiconductor 
material, and the fabrication of multiple concentric quantum rings (up to five) can be 
achieved with high quality and reliability [1]. 

On the other hand, one of the most studied nano-structures, the quantum dot, has 
been appointed as one of the most promising implementations of a single qubit [5l E]. 
There has been a huge amount of work to circumvent the numerous associated problems: 
decoherence [7], the coupling between two qubits (to implement a two-qubit quantum 
gate) in double quantum dots has been studied extensively, in particular how it can be 
tuned using electric fields [8] , magnetic fields [9] , or the effect of the confinement of the 
double quantum dot in a quantum wire [in], etc. 

When dealing with applied magnetic fields, most theoretical studies on quantum 
dots focus on strongly localized states, many times achieved using always-bounding 
potentials, as the three dimensional harmonic potential or impenetrable walls. The 
presence of a constant magnetic field, anyway, is equivalent to a two dimensional 
harmonic potential in the plane orthogonal to the field direction. Conversely, there 
are far less examples of studies considering finite potentials, in particular those whose 
shapes or features allow the presence of resonance states. There is a number of 
reason to study finite potentials, from charge transport situations to the effect of 
resonance states in luminescent quantum dots. The work by Bylicki and Jaskolski 
pTj analyzed the binding of shape resonances through the application of an external 
magnetic field. They considered an one-electron spherical quantum dot-quantum well 
structure (QDQW). In this work we consider a closely related problem with the purpose 
of a better understanding of the transition from a resonance state to a bounded one. 

As has been said above, a constant magnetic field induces a two-dimensional 
harmonic potential that precludes the appearance of resonances in the plane orthogonal 
to the field, i.e. the loss of particles can take place only in the direction of the field, 
closing a number of decaying channels that would be available in absence of the field. In 
this sense, it renders almost irrelevant what particular shape has the binding potential, 
as long as it allows the appearance of resonances in the direction of the field. So, to 
avoid unnecessary complications, we consider a quantum dot model with cylindrical 
symmetry whose axis has the same direction that the magnetic field, besides we use the 
Effective Mass Approximation (EMA). As it is well known in this approximation the 
many-body interactions of the electron trapped in the quantum dot are reduced to a 
(simple) bounding potential, and all the parameters of the Hamiltonian, mass, dielectric 
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constants, and so on are taken as equal to the the bulk parameters. 

Despite its apparent simplicity, the calculation of the resonance states of one 
electron trapped in a bounding potential with an external magnetic field is far from 
trivial, in particular for small enough strengths of the field since the Landau levels are 
bunched when the strength of the field goes to zero. To obtain the energy and width 
of the resonance states we employ the complex scaling method of Moiseyev [12] that, 
together with a variational £^ approximation for the wave function has been used to 
analyze the bound and resonance states of two electron quantum dots [I3]. Moreover, as 
has been shown in [13], the fidelity of the variational eigenstates is a good tool to detect 
the resonance states. In this work we use the fidelity to study the binding of a resonance 
and show that it is signaled by a sharp change in the behavior of the fidelity. This feature 
is consistent with the behavior observed in the fidelity when the system experiments a 
quantum phase transition (in many-body models), is near the ionization threshold or 
to a resonance [13]. Recently, the application of concepts from Quantum Information 
Theory, as the fidelity or entanglement, has been very fruitful to analyze bound and 
resonance states in few body systems as two-electron quantum dots [131 HH [IS [IS], or 
two-electron He-like systems [ITJ [18] . 

The paper is organized as follows, in the quantum dot model and the variational 
approximation that provides approximate eigenvalues and eigenfunctions is presented 
Section [2], the analysis of the binding process is presented in Section [3l and finally we 
summarize and discuss our results in Section HI Some rather lengthly and technical 
results, mostly matrix elements are deferred to the Appendices 

2. Model 

The bounding potential of the quantum dot is given by a piecewise function 

f p<ap, f < l^ol < ^ 

V{v) = I -^2, P < ap, \z\ < f (1) 
[ 0, p>a„ \z\>^ 

i.e. the potential is a cylindrical well aligned with the z-axis, with two potential steps 
at the top and the bottom of the cylinder, where the radius of the cylinder is Cp, its 
height a^, the depth of the potential well, — V2, and the height of the potential steps, Vi. 
This ensures that, for properly chosen constants az^h^^ap^Vi and V2 the one electron 
problem has resonance states without external field. 

Even with the introduction of the external field the angular momentum is a 
conserved quantity, and its eigenvalues good quantum numbers. So, we focus on states 
with zero angular momentum L^. In this case, the Hamiltonian reads as 

-T,icph (4) ^ ^' = ^> = ' 

where V{p, z) is given by Equation [H /i is the effective mass of the electron in the 
semiconductor material. The Equation is written in atomic units. If a magnetic field 
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B = Bz is applied, there is a new term that must be incorporated to the LHS of 
Equation [2] 

i/ = i/„ + |/, (3) 

The discrete spectrum and the resonance states of the model given by Eq. ([3]) can 
be obtained approximately using variational functions [19], [20]. So, if \E'j are the 
exact eigenfunctions of the Hamiltonian, we look for variational approximations 

M 

=^ = E ' = (^^'^)^ ; J = 1, ■ ■ ■ , M . (4) 

1=1 

where the $i must be chosen adequately and M is the basis set size. 

Since we are interested in null angular momentum eigenfunctions, and taking into 
account the symmetries of the Hamiltonian |3l we choose as basis functions 

^i{p,z) =ipn{vp)M'^z)^ (5) 



where 



and 



Mvp) = ^=r/e-'"'/2£W (r/p) , (6) 
Vn + 1 



'/'.M = ^^e-^'^'/'/:r(H^I)), (7) 

7] and u are the variational parameters. The matrix elements of the kinetic energy, the 
bounding potential and the magnetic field term are given in the Appendix. With all 
these matrix elements we get a variational eigenvalue problem 

Hc(^) = E]c'^^\ (8) 

where the entries of the matrix H are given by 

ilk,n,l,m = {tpk(pn\H\lpl(f>m) (9) 

The ample range of materials and structures available to design self-assembled 
quantum dots precludes the possibility of a very general analysis, however, in part for 
comparison reasons, we use similar parameters than those used by Bylicki in [21] and 
Bilicky and Jaskolski in [TT] to model a quantum dot quantum well structure made of 
gallium arsenide (GaAs) composites, i.e. ap = = 7nm, = 2.5nm, V\ = 0.37eV, 
Vi = 0.108844ey, and /i = 0.041mo. The potential well depth is slightly different from 
the value used by Bylicky, but should be chosen so owed to the cylindrical symmetry 
employed in this work, conversely to the spherical one used by Bylicki. As we will show, 
this set of parameter is consistent with a resonance energy around 20meV with B = 0. 

Figure [1] shows the variational eigenvalues obtained for V2 = 0.108844, as functions 
of the magnetic field strength B. As can be appreciated, the eigenvalues are grouped in 
sets that are bounded between Landau levels given by Ell^M^h) = — (|M| + n + 1/2), 
where M is the azimuthal angular momentum ( in our case M = 0), and n = 0,l,2,.... 
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Figure 1. (Color on-line) The four panels show details of the variational spectrum 
and the Landau levels, a) The variational spectrum obtained with a basis set size 
= 2500. Each black curve corresponds to a single eigenvalue and the red curves 
correspond to the exact Landau levels, b) This panel shows the limit between the zone 
limited by the two first Landau levels and the zone above the second Landau level. 
The only red curve corresponds to the second Landau level. The figure shows how the 
levels that collision with the second LL have several avoided crossings, c) The difference 
between the lowest variational eigenvalue and the first LL, when the lowest eigenvalue 
crosses below the first LL the eigenstate becomes localized, d) The eigenvalues above 
the first LL. The figure shows clearly how the eigenvalues are, basically, parallels to 
the LL and accumulate above it. 



A pair of Landau levels delimits a region of the (S, E) plane, in this zone the eigenvalues 
are parallel to the lower Landau level, except for the appearance of avoided crossings. 

As can be seen very clearly from Figure [1], the variational eigenvalues have the 
tendency to accumulate above the Landau levels. This feature is similar to the eigenvalue 
accumulation observed above the continuum threshold [I3]. When a magnetic field is 
applied, each Landau level works as the bottom of a continuum leading to the eigenvalue 
accumulation observed in Figure [H Another salient feature of the spectrum appears for 
large enough magnetic fields: one isolated eigenvalue with lower energy than the lowest 
Landau level. From physical point of view, the origin of this state can be understood as 
follows: for intermediate values of the magnetic fields the wave function of the electron 
looks like an harmonic oscillator wave function in the (p, 0) plane and as a free particle 
in the z direction. The spatial extent of the wave function on the plane is roughly 
equivalent to the radius of a Landau level, that is larger than the radius (in the plane) 
of the quantum dot. When the strength of the magnetic field is increased the radius 
of the lowest Landau level becomes smaller and smaller reaching, at some point, a size 
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similar to the radius of the quantum dot, at this point the state becomes locahzed. As 
we will show, the mechanism of localization can be quantified and strongly influences 
the behaviour of the resonance states whose energy lies near the localization point. The 
spectrum in Figured] was obtained using a basis set size N = 2500, all the other results 
shown in what follows were obtained using = 3200. 

3. Detecting the resonance states using complex exterior scaling 

The calculation of the energy and width of resonance states offers a number of challenges 
that repeatedly leads to the formulation of new methods. Among the most widely used 
methods can be mentioned the complex scaling (or complex dilation) method [12], the 
complex absorbing potential method [22] and the density of states method [23]. Each 
one of these has its advantages and drawbacks. Because we are dealing with a piecewise 
potential we resort to the exterior complex scaling of Moiseyev [12]. This method is 
particularly appealing in our case since for 5 7^ 0, the complex scaling should be applied 
only to the z variable because this is the only direction available to the electron to get 
away from the bounding potential. The exterior of the method refers to the exterior of 
the region where the potential is not trivial so, for B ^ the exterior complex scaling 
asks that 



where z is the coordinate to be complex scaled, and 6 is the rotation angle. The complex 
scaling turns the non-normalizable scattering states of the continuum and the resonance 
states into normalizable states (with a different norm) which can be analyzed with the 
full variational approach usually employed in Hermitian problems. 

Another advantage of the complex scaling method lies in the fact that its 
implementation depends on the availability of matrix elements that are similar to the 
ones calculated to obtain the entries of iik^n,i,m- 

Figure [2] shows the behaviour of the energy and width of the lowest resonance as 
a function of the magnetic field strength. The binding of the resonance can be clearly 
observed, i.e. the width of the resonance drops to zero when the magnetic field strength 
increases its value. The exterior complex scaling methods works nicely for intermediate 
values of the magnetic field. Anyway, if i? < 5T, or i? > 16T, the convergence of 
the method is, at least, questionable. For small fields the method struggles to provide a 
reliable value for the resonance eigenvalue because, as can be appreciated from Figure El 
ioT B ^ 5 the resonance "enters" in the continuum above the second Landau level (LI). 
In this region the approximate eigenvalues that enter from the region above the first 
Landau level "collide" with the eigenvalues that lie between the first and second LL. The 
successive collisions, and the corresponding avoided crossings, can be well appreciated 
in Figure [Tb). In this sense the resonance disappears when there are not open channels 
to decay, bellow the first LL there is only a localized state, so there is no way for the 
resonance to avoid the binding. 
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Figure 2. The real and imaginary parts of the resonance eigenvalue E{9) vs the 
magnetic field strength, a) The solid square (blue) dots correspond to the energy 
of the resonance and the two solid lines correspond to the two lowest Landau level 
energies. The value for B = was obtained performing the complex rotation in the 
coordinates p and z. b) The width of the resonance vs the magnetic field. As the 
magnetic field increases its value the width of the resonance goes to zero, signaling its 
binding. 



The scenario between LL's can be better appreciated in Figure [31 The Figure shows 
the complex spectrum obtained when the complex rotation is performed accordingly 
with Equation [in] for different values of 6. In this Figure it is clear why the method is 
termed complex rotation, the continuum part of the spectrum now lies in the complex 
plane over straight lines, the angle between the lines and the real axis equal to 26. When 
B ^ 0, the data for different ^'s form "hand-fans" of data, i.e. sets of straight lines with 
a common origin in the real axis. Each one of these sets can be associated to a single 
LL, the leftmost hand-fan corresponds to first LL, the middle one to the second LL, and 
the rightmost to the third LL. Again, this structure can be attributed to role that each 
LL plays as the bottom of a continuum of levels. Effectively, the distance between the 
origins of each set of lines, over the real axis, is exactly equal to the distance between the 
LL at the magnetic field strength considered. As the magnetic field strength decreases, 
the data of the different fans overlaps making very difficult to distinguish the eigenvalues 
corresponding to the resonance. 

Figure [3b) shows the trajectory in the complex plane of the eigenvalue associated 
to the resonance state. It is well known [12] that, to look for the best value obtainable 
from the method, the complex rotation must be performed for different values of 6 and 
the best approximation corresponds to the stationary points of the trajectory. Actually, 
this is the recipe that we followed to obtain the data shown in Figure [2^) and b). It 
is possible to modify the method to obtain approximate resonance energies and widths 
for B < 5T, we will be back at this point latter on. 
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Figure 3. (Color on-line) a) The complex rotation eigenvalues vs the magnetic field. 
From top to bottom, the eigenvalues were obtained for 9 = 0.1 (black • ), = 0.11 
(red ■), = 0.12 (green ♦), up to 6* = 0.2. The three fan-like sets of data are related 
to the first three LL's. Is is clear that the the leftmost fan and the central one start 
to overlap around 27 meV. For B < 5T the overlap grows larger making more difficult 
to look for the resonance data, b) The 9 trajectory of the resonance eigenvalue. The 
symbols correspond to the angles shown in panel a). This eigenvalues can be observed 
in panel a), they lie near the bump that presents the leftmost fan around 21 meV 



3.1. Localization 

The binding of the resonance for large enough magnetic field is a consequence of 
the strong localization experimented by the wave function. As a matter of fact, the 
localization is so sudden and clear that allows to follow the resonance through the 
spectrum. To quantify the localization of an approximate eigenfuntion \E'J, we calculate 
the probability that the electron is localized in the potential well, 

Pj= r r \'^ap^^)? dpdz, (11) 

J-^^ Jo 

where Pj is the probability attributable to the localization of 

Figure H^) shows the probability Pi as a function of the magnetic field strength 
for the lowest variational eigenvalue, while Figure |4)d) shows the probability Pj for 
i = 2, 3, . . . , 8. The curves are shown in different panels because of their respective 
scale. 

It is worth to mention that despite the variational eigenvalues near the localization 
do not show any sudden changes, except for the lowest one when it crosses the first 
LL, the probability of localization into the well is appreciable for a large number of 
eigenvalues. As we will shown later, the localization probability allows to follow the 
presence of the resonance from the localization threshold up to i? ~ 5T. 




B[T] B[T] 

Figure 4. The localization probability vs the magnetic field strength. The left panel 
shows Pi, i.e. the probability corresponding to the lowest variational eigenvalue. 
For small enough magnetic fields the variational eigenstate is extended over a region 
much larger than the characteristic lengths of the quantum dots. When B > 17T 
the eigenstate becomes localized, as shown by the sudden increase of Pi. Once 
the lowest eigenstate becomes localized it remains so for even larger field strengths. 
Right panel: From top to bottom, P2 (solid line), P3 (dotted-line), . . ., Pg. The 
localization probabilities show a peak where the corresponding eigenvalue crosses near 
the resonance, after and before the eigenstate is extended. 



3.2. Fidelity 

For a given quantum state ip, that depends on a parameter A, a measure of how much 
it changes when the parameter is varied is given by the fidehty, J-", which is defined as 

-^aa(A) = m\ - AA),7A(A + AA))|' , (12) 

where AA is a small variation of the parameter. The fidelity has been extensively 
used to characterize the analytical properties of quantum states near a quantum phase 
transition in quantum spin chains models [23], quantum phases of matter [25], or bound 
and resonance states in atomic or quantum dot models [13]. In two-electron quantum 
dots, it has been shown that the fidelity of the approximate variational eigenstates detect 
the resonance states and allow to calculate approximately its energy [13]. To achieve 
this, the fidelity of many eigenstates should be calculated as a function of the external 
parameter that drives the system from bound to unbounded states. In this work, this 
parameter is the magnetic field strength. 

From its definition. Equation [121 it is clear that most of the time J-aa ~ 1, except 
for very special cases. On the other hand the fidelity should drop to zero if the system 
experiments a rather sudden change [21]. For this reason, often it is convenient to study 
the function 

^„(A) = 1 - - AA), vl;;;(A + AA))|^ . (13) 

Figure |5] shows the behaviour of the Qn function for n = 1, 2, . . . , 6, as function 
of the magnetic field strength. Each curve has a more or less well defined peak for 
a given value of the magnetic field, -B„. The values En{Bn), i.e. the value of the 
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Figure 5. The function C/„ vs the magnetic field. The figure shows the function for 
the first six variational eigenvalues, Qi , G2, ■ ■ ■, Ge, solid ( ), dotted ( ), short- 
dashed ( ), long-dashed ( ), dot-dashed ( — • — ) and double dot-dashed 

( — • • — ) lines respectively. The peak where each curve attains its maximum value is 
clearly appreciable. i?„ is given by the abscissa of the peak (see the text). 



variational eigenvalues at their respective peaks of the function Qn, give a very good 
approximation for the resonance energy at the points As has been pointed out [13] 
this does not allow to obtain the resonance energy for every magnetic field strength 
value, nevertheless the method provides another tool to analyze resonance states. The 
fidelity method works best when the width of the resonance is not too wide, so it is to 
be expected that it will be more precise near the localization point. 

Before presenting the results of the localization and fidelity methods with respect to 
the resonant states, we want to stress the relationship between both quantities. Let us 
call Bp{ap) the critical value of the magnetic field such that the localization probability 
of the first variational eigenvalue becomes noticeable for a given quantum dot radius 
Op, and Bjr[ap) the magnetic field value such that the function Qi attains its maximum 
value as a function of B. Figure |6] shows the values of both quantities, Bp{ap) and 
Bjr[ap), for several values of the quantum dot radius ap. The agreement between both 
critical quantities is striking. 

At this point, we can summarize our results in Figure [71 The Figure shows the 
energy of the resonant state as a function of -B, as it is obtained from the complex 
rotation, fidelity and localization methods. 

From Figure [7] is rather clear that the localization methods is able to follow the 
resonance from the localization point until the second Landau level. The same can not 
be said about the fidelity results that, anyway, it is known that are valid for resonance 
with very small widths. 

Above the second LL the complex rotation method should be complemented with 
some ad hoc assumptions to identify the points that belong to the ^-trajectory. The main 
assumption is rather physical: the complex rotated eigenvalues can be approximated by 
Ell + tan2^ii^„, i.e. they form a straight line in the complex plane that form an angle 
of 26 with the x-axis. En is a variational eigenvalue. This approximation works nicely 
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Figure 6. (Color on-line) The critical fields Bp{ap) and Bjr[ap) vs the quantum dot 
radius. The localization probability critical field (red squared dots, ■) and the fidelity 
critical field (solid black dots, • ) data is shown for several quantum dot radius. The 
(blue) dashed curve correspond to the radius of the lowest LL as a function of the 
magnetic field. It is clear that for small quantum dot radius the localization takes 
place when the cylindrical wave function enters into the cylindrical QD. For larger 
QD's radius the localization is dominated by the quantum well potential and not by 
the magnetic field. 



except near a resonance. The second assumption is similar: a given complex eigenvalue 
E^{6) should not change too much if 6 is changed, say 



where SO is small enough. Applying these assumptions and discarding complex 
eigenvalues that jump from a hand-fan to another (see Figure EK)) it is feasible to 
obtain stabilized values for Eres above the second LL. 

4. Discussion and Conclusions 

The external magnetic field precludes the escape of the electron on the (p, (p) plane, 
for this reason the quantum dot bounding potential consists of a potential well and a 
potential step, to ensure the presence of resonance states. In this sense, the problem 
has several characteristic lengths, the quantum dot radius, the Landau levels radii, the 
length of the quantum dot along the z axis, and so on. The two lengths that come into 
play varying the magnetic field strength, for the set of parameters considered in this 
work, are the quantum dot radius and the lowest LL radius. The potential well depth 
and the other parameters of the QD were chosen to ensure that for B = there were 
no bound states. 

The localization probability is able to track down the resonance states from the 
localization point up to the second LL. Above the second LL the eigenvalues have too 
many avoided crossings, rendering the method useless. The many avoided crossings 
also prevent the use of the fidelity to detect the resonance above the second LL. The 
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Figure 7. The resonance energy Eres calculated using the three methods depicted in 
the text, complex exterior scaling, localization probability, fidelity, the corresponding 
values are shown using blue square dots (■), black circles (O ), and magenta diamonds 
(<)). The energy of the resonance in the region above the second LL can be obtained 
using a modified version of the complex rotation method and is shown using blue open 
squares dots (□) . 



large width of the resonance between the first and second LL seems to severely limit 
the ability of the fidelity to detect the resonance in this region. In two-electron dots 
the range of utility of the fidelity method could be related to the expectation value of 
the Coulomb repulsion between the electrons |T3j, right now we do not have a simple 
physical argument to explain the behaviour observed in this work. 

Finally, both critical fields, Bp{ap) and Bjr{ap), show two clearly distinguishable 
regimes, for small quantum dot radius their values are given by the lowest Landau 
level radius, so both are proportional to l/a^. For large enough quantum dot radius 
both critical fields show a different behaviour and, apparently, both are proportional 
to (a — PttpY, where a y /3 are constants. The extent of the small field regime can 
be tuned changing the parameters of the bounding potential [H extending or reducing 
it. The critical behaviour of the eigenvalues in the transition region between the large 
field-small quantum dot radius and the small field-large quantum dot radius will be 
analyzed elsewhere. 
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Appendix A. Some matrix elements 

It is convenient to separate the kinetic energy matrix elements in two contributions, one 
corresponding to the radial coordinate. 
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(V'n \Tp\ iIJs) = (i^n 



V 



and the other corresponding to the z coordinate, 

1 



where 



while 



— 



2^^{n + l){s + l) V4 2 ^ 



Ti = (n + 



n s—1 



p=0 9=0 
n— 1 s 

p=0 13=0 
n— 1 s—1 



p=0 9=0 



(n- 1 - p)\{2 + p)\p\{s - q)\{l + q)\qV 

[-l)p+i[n + + l)!(p + g + 1)! 
(n - p) ! (1 + p) !p! (s - 1 - g) ! (2 + g) !g! ' 

(-l)P+g(n+l)!(g + l)!(p + g + l)! 
[n-l - p)\{2 + p)\p\{s - I - q)\{2 + q)\q\ 



T2. 



Iz — ^t,r^ 

t r-1 

EE 

d=0 /=0 
i-1 r 

EE 

d=o /=o 

t-l r-1 

EE 

d=0 f=0 



-iy+H\ r\{d + f)\ 



3z 



4z 



^t-d)\id\yir-l-f)\il + f)lfV 

(-iy+rt\ r\{d + f)\ 
(t-l-rf)!(l + rf)!(r-/)!(/!)2' 

(_l)t+r^! r!(rf + /)! 
(t_l_rf)!(l + c?)W!(r-l-/)!(l + /)!/!• 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 
(A.7) 



The matrix element of the bounding potential can be factorized owed to its piecewise 
character and using that there is a potential barrier on the z direction. We get that 



{lpn4>t \V{p, Z)\ij,(t)r) = -V2 Iv2 + Vl hl^-n 



the barrier term is, obviously, proportional to Vi, while the term proportional to V2 
corresponds to the matrix element of the potential well. The matrix elements Iv2 and 
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lyi are given by 

-l)P+''(n + l)!(g + l)! Iv2i 
V(n + l)(s + l) ^ ^ - P)K1 + - + 

(t-rf)!(rf!)2(r -/)!(/!) 



d=o /=o 



and 



{p + q + l- k) 



fc=0 



fe=0 



{d+f-k)\ 



The matrix element of the magnetic field term reads as 

{ipn \Hc\ ^s) 



Wv^(n + + 1) 



n s 
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